Anatomy of a Merger: A Numerical Model of A754 
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ABSTRACT 

A754 is a well-observed cluster of galaxies (z=0.054) which exhibits a variety of 
morphological peculiarities. These include a bar of X-ray emission that is offset signifi- 
cantly from the galaxy distribution, an elongated X-ray surface brightness distribution 
extending between two distinct peaks in the galaxy distribution, and an extremely non- 
isothermal and asymmetric intracluster medium (ICM) temperature morphology. Using 
these observational constraints, we present a numerical Hydro/N-body model of A754 
in which two clusters (2.5:1 mass ratio) have merged nearly in the plane of the sky less 
than 0.5 Gyrs ago with an impact parameter of ~120 kpc, and an impact velocity of 
~2500 km s~^ (roughly the escape velocity of the primary cluster). Our models allow 
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us to identify the origin of A754's peculiar X-ray and temperature morphologies, the 
underlying hydrodynamical processes that shape them, and their future evolution. We 
make detailed predictions for future high resolution X-ray spectroscopic observations 
(e.g. ASTRO-E). We discuss general properties of our models which will be character- 
istic of off-axis mergers. In particular, we find significant non-thermal pressure support 
within the central region which could bias cluster mass estimates. We find significant 
angular momentum imparted on the gas distribution in the cluster. We find that mixing 
of the subcluster gas components is an inefficient process, particularly at large radii. Fi- 
nally, we find that subsequent dark matter core passages result in an extended relaxation 
timescale. 

Subject headings: Galaxies: clusters: individual (A754) - galaxies: intergalactic medium 
- hydrodynamics - methods: numerical - X-rays: galaxies 
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1. INTRODUCTION 

In recent years, clusters of galaxies have become 
essential cosmological probes. Both their current dy- 
namical state and internal structure provide impor- 
tant clues to conditions in the early universe. Over 
the years, many researchers (e.g., Geller & Beers 1982; 
Dressier & Shectman 1988; West & Botliun 1990; 
Jones & Forman 1991; Bird 1994; Davis 1994; Mohr 
et al. 1995, among others) using both X-ray and op- 
tical databases, have attempted to understand the 
dynamical state of clusters through analysis of their 
projected substructure (i.e., dumpiness and asymme- 
tries in the distribution of the intracluster medium 
(ICM) and galaxies). Unfortunately, these studies 
are inherently limited because they represent only 
a 2-dimensional projection of the true 3-dimensional 
cluster shape, and moreover, there is no information 
regarding the gas kinematics. As such, the current 
data set cannot reveal the underlying physical pro- 
cesses which produce the observed substructure nor 
can they reveal how the substructure will evolve in 
time. For these answers, one must turn to numerical 
simulations (e.g., Evrard 1990; Cen 1992; Roettiger 
et al. 1993; Schindler & Miiller 1993; Pearce, Thomas 
& Couchman 1994; Navarro & White 1994; Mohr et 
al. 1995, among others). In this paper, we apply nu- 
merical hydrodynamical/N-body simulations to the 
analysis of the extensive observational database col- 
lected on the galaxy cluster Abell 754. 

A754, at z=0.054, is one of the most extensively ob- 
served clusters in both the optical and X-ray bands. 
Much attention has been drawn to this particular 
cluster by its peculiar X-ray morphology (Fabricant 
et al. 1986, F86; Henry & Briel 1995, HB95). Far from 
being a relaxed spherical distribution centered on the 
galaxies, A754 shows an elongated and barred X-ray 
morphology significantly offset from a bimodal galaxy 
distribution. Even more recently, A754 has been rec- 
ognized for its significant non-isothermality (IIB95; 
Henriksen & Markevitch 1996, hereafter IIM96). To- 
gether, these properties make A754 an ideal labora- 
tory for the study of non-equilibrium systems. 

The purpose of this paper is to assemble the ob- 
servational database on A754 and create a plausible, 
though not necessarily unique, self-consistent, numer- 
ical model that is in general agreement with the obser- 
vational data. The model is then used to explore the 
underlying hydrodynamical processes that produced 
the observed morphological features, with particular 



emphasis on the X-ray surface brightness and temper- 
ature distribution. We then follow the development 
of these features as the merger is allowed to evolve 
thus giving insight into the more general evolution of 
a slightly offaxis merger. 

We perform this study using a hydrodynamical/N- 
body code based on the piecewise parabolic method 
(PPM) and a particle-mesh algorithm (PM). The 
gaseous intracluster medium (ICM, as represented by 
PPM) is allowed to evolve self-consistently (including 
self-gravity) within a changing gravitational poten- 
tial defined largely by the dark matter distribution as 
represented by the coUisionless N-body particles. In 
this study, we employ idealized initial conditions in a 
manner similar to previous studies by Roettiger et 
al. (1993, 1995, 1996, 1997a), Pearce et al. (1994), 
and Ricker (1997). The Pearce et al. simulations are 
based on a Smoothed Particle Hydrodynamics code 
(SPH, Lucy 1977), while Ricker (1997) explores more 
extremely offaxis mergers. The idealized initial con- 
ditions used here allow us to explicitly control the 
merger parameters making possible a systematic sur- 
vey of parameter space. It also allows for an efficient 
use of the computational volume providing improved 
resolution over simulations of largescale structure for- 
mation. On the other hand, our study is limited by 
the neglect of largescale tidal forces, cosmological in- 
fall, and uncertainty in the detailed structure of clus- 
ters. The cosmological infall of baryonic matter is 
not likely to affect the inner regions of the cluster 
which are of interest here. Of potentially greater im- 
portance, these simulations do not include radiative 
cooling which may influence the cluster's thermal evo- 
lution, particularly in the high gas density core. 

In §|2] we review the recent optical and X-ray ob- 
servations of A754. Section ^ describes the numerical 
method, the code, grid configuration, etc. In § ^ we 
discuss the degrees of freedom in the model and the 
manner in which the observations are used to con- 
strain the model parameters. Here, we also describe 
the initial conditions and summarize the merger pa- 
rameters. In we make a quantitative comparison 
of our model with the observational data while dis- 
cussing the current dynamical state of A754. Section 
describes the future evolution of the system. We 
summarize our results in S|7] We define h =Ho /lOO for 
consistency with previous observational results, but 
we choose Ho =65 km s~^ Mpc~^ for the scaling of 
our model. 
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2. THE OBSERVATIONAL DATA 

2.1. Optical Data 

Many studies have been performed on tlie distribu- 
tion of galaxies in A754. Of these, many have found 
significant substructure (e.g., Geller & Beers 1982; 
F86; Escalera & Mazure 1992 (EM92); Bird 1994; and 
Zabludoff & Zaritsky 1995 (ZZ95)) while others have 
not (e.g.. Dressier & Shectman 1988, West & Bothun 
1990). Of those studies that do find substructure, two 
primary galaxy concentrations are consistently iden- 
tified. Their spatial location relative to the X-ray 
surface brightness distribution is indicated in Figure 
H (see also Plate L8 in ZZ95). Hereafter, we refer to 
these as the SE and NW concentrations. Table 1 is 
a summary of the individual subcluster properties as 
determined by several researchers. 

Examination of Table 1 reveals a general consis- 
tency between the various analyses. The largest dis- 
crepancy resides in the velocity dispersion of the NW 
population determined by EM92. Their value differs 
from the other two by greater than 2ct and is less than 
the dispersion calculated for the SE group. In each 
of the other two studies, the NW group has a larger 
velocity dispersion than the SE group. It is also in- 
teresting to note that in each study more galaxies are 
found to be associated with the SE subcluster than 
with the NW subcluster. 

The various estimates of A754's virial mass are 
consistent although uncertainties are large. Biviano 
et al. (1993) find 4.07±|;0| x lO" Mq and 6.03±|67 
X 10^'* Mq within 0.75/1"^ and l.bh'^ Mpc, respec- 
tively. Escalera et al. (1994) find 10.48±3.27 x 10^^ 
Mq corresponding to a virial radius of 2.21±0.28/i~^ 
Mpc. Only EM92 have attempted to determine 
masses for the individual concentrations. They find 
2.9 X 10" Mq and 1.3 x lO" Mq within -0.32 Mpc 
for the SE and NW concentrations, respectively. 

To summarize, the data are consistent with two 
clusters at nearly identical redshifts having a pro- 
jected separation of O.TSh^^ Mpc (ZZ95). Their ve- 
locity dispersions are of order 900 km s^^ with the 
ratio of their dispersions consistent with unity. The 
total mass of the system within ^2h^^ Mpc, based on 
the galaxy dynamics, is likely greater than 10^^ Mq. 

2.2. X-ray Data 

The ROSAT X-ray surface brightness image gen- 
erated by IIB95 reveals an unusually complex mor- 



phology. We have reproduced this data in Figure 
|l| (solid contour). There are several features worth 
noting. First, the general elongation which extends 
parallel to the line connecting the galaxy concentra- 
tions (§2.1.). Second, the peak of the X-ray distri- 
bution is seen to be a well-defined bar which extends 
up to ^A' to the north of the SE galaxy concentra- 
tion and is nearly perpendicular to the line connect- 
ing the galaxy concentrations. In a wavelet analysis 
of the X-ray surface brightness distribution (Slezak, 
Durret, & Gerbal 1994, Figure 16 therein.), the peak 
of the distribution is associated with an angled feature 
which includes the bar and a less well-defined exten- 
sion toward the NW galaxy concentration. The total 
effect is that the eastern edge of the bar is sharper 
than the western edge. Finally, it is important to 
note that neither of the major galaxy concentrations 
are spatially coincident with the peak X-ray emission 
(ZZ95). The dominant galaxy in A754, identified as 
a cD galaxy (Dressier 1980), is located in the NW 
subcluster, 13' from the X-ray peak. 

The mean temperature of the X-ray emitting gas 
in A754 is 9 keV (IIM96). Using the global veloc- 
ity dispersion of 750 km s~^ (Mazure et al. 1996), 
(3 = ^mh(Jy/kT is found to be 0.4 which is signif- 
icantly discrepant from the best-fit mean value of 
0.94±0.08 found by Lubin & Bahcall (1993) for a large 
sample of nearby clusters. In any case, the mean tem- 
perature may be misleading in that A754 is decidedly 
non-isothermal. Henriksen (1993), using HEAO-l A2 
data, identified two ill-defined temperature compo- 
nents within A754. IIB95, using ROSAT data, were 
the first to produce a detailed temperature map of 
A754. They found a relatively cool region near the 
X-ray peak at ~5-7 keV. To the north and south they 
identify arcs of gas greater than 9 and 10.8 keV, re- 
spectively. To the west, they identify an arc with 
temperature greater than 13.8 keV. 

A similar attempt to map the temperature distri- 
bution was made by HM96 using ASCA data. Al- 
though differing in detail, there is considerable con- 
sistency between the HB95 and IIM96 maps. Only 
one of the nine central regions (as defined by HM96, 
Figure 1, therein) differed by a statistically signifi- 
cant amount from the corresponding region in IIB95. 
Even in the low surface brightness outer regions of 
the cluster where the detailed energy dependence of 
the ASCA point spread function and the background 
subtraction are both important, there is considerable 
agreement between the two maps. It is likely that 
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much of the observed discrepancy results from the 
different energy bands employed. The low energy 
ROSAT band (0.1-2.5 keV) is less able to constrain 
the temperatures of the extremely hot gas observed 
by ASCA which is sensitive to the 0.5-10 keV band. 

Finally, for completeness, we comment on the ev- 
idence for a cooling flow in A754. Edge, Stewart & 
Fabian (1992) report evidence for a moderate cool- 
ing flow, ~24 M0/yr. However, there is no evidence 
for the optical filaments often associated with cooling 
flows (Heckman et al. 1989). Both HB95 and HM96 
report relatively low temperatures (~6 keV) near the 
X-ray maximum. It as been suggested that this is cool 
gas having been stripped from one of the subclusters 
(HM96). It is also possible that this is a remnant of 
a cooling flow that was disrupted by the merger. No 
attempt was made to model a premerger cooling flow 
although we would expect one to have been disrupted 
during the merger (Gomez et al. 1997). If there had 
been a long-established cooling flow in A754 before 
the merger, it is quite likely that remnant cool gas 
would remain at this early epoch of the merger ow- 
ing to the inefficient mixing of cluster gas components 
(See §|3]). 



3. NUMERICAL METHOD 

We compute the dynamical evolution of merging 
clusters of galaxies using a hybrid hydrodynamical/N- 
body code in which the hydrodynamical component is 
CMHOG written by one of us, J. M. Stone. CMHOG 
solves the fluid equations using an implementation 
of the piecewise-parabolic method (PPM, Colella & 
Woodward 1984) in its Lagrangian remap formu- 
lation: its fidelity has been demonstrated through 
numerous applications in interstellar gas dynamics, 
e.g., the interaction of strong shocks with dense clouds 
(Stone & Norman 1992). The coUisionless dark mat- 
ter is evolved using an N-body code based on a stan- 
dard particle-mesh algorithm (PM, Hockney & East- 
wood 1988). The particles are evolved on the same 
grid as the gas using the same time step. The time 
step is determined by applying the Courant condition 
simultaneously to both the dark matter and the hy- 
drodynamics. The only interaction between the col- 
lisionless particles and the gas is gravitational. Since 
we are modeling an isolated region, the boundary 
conditions for Poisson's equation are determined by 
a multipolc expansion of the mass distribution con- 
tained within the grid (Jackson 1975). Particles that 



leave the grid are lost to the simulation. Typically, 
less than a few percent of the particles leave the 
grid. Previously, this code has been used to examine 
systematic errors in the measurement of the Hubble 
constant using Sunyaev-Zeldovich effect in merging 
clusters of galaxies (Roettiger, Stone, & Mushotzky 
1997). 

The simulation is fully three-dimensional. Two dif- 
ferent computational grids were used. A survey of 
parameter space was performed on the MasPar-2 at 
Goddard Space Fhght Center (GSFC) using a fixed 
and rectangular grid (256 x 128 x 128 zones) hav- 
ing linear dimensions of 12.8 x 8.3 x 8.3 Mpc. The 
merger axis coincides with the grid's major axis along 
which resolution is uniform and scales to 50 kpc or 
^5 zones per primary cluster core radius. Resolution 
along the grid's minor axis is uniform within the cen- 
tral 64 zones and ratioed in the 32 zones on either 
side. That is, the resolution is a uniform 50 kpc ex- 
tending 1.6 Mpc (32 zones) on either side of the major 
axis. Beyond 1.6 Mpc, the zone dimensions increase 
by ~3% from one zone to the next out to the edge 
of the grid. We use outflow boundary conditions for 
the hydrodynamical evolution. Once we settled on a 
set of initial conditions, we ran a very high resolution 
simulation on the CM-5 at the National Center for 
Supercomputing Applications (NCSA) using a fixed 
rectangular grid with dimensions of 512 x 256 x 256 
giving a resolution of 25 kpc in the central 3.2 Mpc 
(128 zones), or nearly 10 zones per core radius. 

4. MODEL CONSTRAINTS 
4.1. The Data and its Limitations 

We use the above observations to constrain our nu- 
merical model. However, there are serious limitations 
to this process. First, the observational uncertain- 
ties can be large. Second, we believe, as did ZZ95, 
that this is a major merger in the very early stages 
after core passage. If this is the case, then the struc- 
ture, temperature and even the mass of the initial 
systems can be severely obscured by the current non- 
equilibrium conditions (Roettiger et al. 1996). Fi- 
nally, there are just too many free and not necessarily 
independent parameters to completely constrain the 
model. These include, the total mass of the system, 
the relative masses of the individual subclusters, the 
distribution of mass within the subclusters (p ^ r~", 
core sizes, etc), the fraction of mass contained in 
baryons, the distribution of baryons relative to the 
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dark matter (i.e., the fi parameter, central gas den- 
sity), the angular momentum of the system (impact 
parameter and velocity) , the linear scaling (we assume 
Ho =65 km s~^ Mpc~^), the merger epoch (time since 
closest approach), and the orientation of the merger 
with respect to the observer. For these reasons, it 
is difficult to claim a definitive model. Rather, we 
have produced a model which agrees well enough with 
the observational data that we may consider it to be 
a plausible representation of the gasdynamics within 
A754. 

Our approach is to use the observable properties of 
the general cluster population to construct our initial 
clusters. We then explore essentially two branches 
of parameter space (rela tive central gas density and 
impact parameter, § 4.3. ) at low resolution (50 kpc 
per grid zone), determine a best fit, then run a high 
resolution (25 kpc per grid zone) simulation which is 
discussed in §^ 

4.2. Initial Conditions 

Since it is likely that much of the information about 
the internal structure of the premerger clusters has 
been erased by the merger, we begin with clusters 
that are consistent with observations of presumably 
relaxed systems. Both subclusters were chosen to 
have the same mass distribution. For our clusters, 
we have chosen the lowered isothermal King Model 
described in Binney and Tremaine (1987). The low- 
ered isothermal King Model is a family of mass distri- 
butions characterized by the quantity ■f/'/cr^ which es- 
sentially defines the concentration of matter. As V'/c^ 
increases, the concentration parameter also increases, 
i.e., the core radius (r^,) decreases relative to the tidal 
radius, r^. We have chosen a model with ■0/(7^=12 
in which we have truncated the density distribution 
at 15rc. Near the half mass radius, the total mass 
density follows a power law distribution, p = r~" 
where a 2.6. This model is consistent with mass 
distributions produced by the cosmological N-body 
simulations of Crone et al. (1994) which showed a ~ 
2.4 in a high density universe (f^=l) and a ^ 2.9 in a 
low density universe ($7=0.2). Observationally, galax- 
ies are distributed as a ^ 2.4±0.2 (Bahcall & Lubin 
1994). 

The simulations are conducted in scaled units. 
They are rescaled to meaningful physical dimensions 
after the fact. This gives us some flexibility in defining 
the initial conditions in that dimensionless subcluster 
parameters are more important than the global scal- 



ing. The observational data does give some insight 
into the scaling of global parameters. The various 
mass estimates listed in S2.1. seem to indicate that 



the total mass within 2-3 Mpc is likely greater than 
1.0 xlO^^ Mq. The observed redshift of A754 indi- 
cates that the galaxy concentrations are separated by 
~1.12 Mpc for Ho =65 km s'^ Mpc"!. 

The attempts to separate the two galaxy compo- 
nents can be used to constrain their relative masses. 
There is considerable uncertainty in doing so since at 
this stage of the merger the velocity dispersions may 
not be representative of the individual masses. We 
have found in previous (e.g., Roettiger et al. 1993) 
simulations that subclusters are "heated" relative to 
the primary as they pass through or near its core. 
Consequently, they may appear more massive than 
they actually are. Also, as seen in Table 1, small 
number statistics lead to large uncertainties in the 
velocity dispersions. The ratio of the NW to SE sub- 
cluster velocity dispersions are consistent with unity 
which would indicate that they are of comparable 
mass (assuming heating is not significant). Taken at 
face value, the velocity dispersions in F86 and ZZ95 
indicate that the NW subcluster is somewhat more 
massive than the SE subcluster, while EM92 indicates 
that the SE subcluster is more massive. Looking at 
the galaxy counts in Table 1 and assuming that the 
galaxies trace the dark matter mass and that signifi- 
cant exchange of galaxies has not occurred, all three 
samples in Table 1 indicate that the SE subcluster 
is the more massive subcluster. HB95 also suggests 
that the SE subcluster is the original main cluster 
based on its central location relative to the presum- 
ably undisturbed outer X-ray contours. For these rea- 
sons, we consider mergers where the SE subcluster is 
the more massive cluster by a factor of 2.5. Although 
marginally consistent with the velocity dispersions of 
F86 and ZZ96, it does represent one extreme, and a 
lower mass ratio (approaching 1:1) cannot be ruled 
out. 

4.3. Parameter Space 

Having settled on the relative cluster masses, we 
then explored a range of impact parameters and rel- 
ative central gas densities. We examined impact pa- 
rameters ranging from zero (head-on) to the sum of 
the respective core radii (Table 2). We found that in 
order to generate the combination of compression and 
shear necessary to produce the X-ray core morphol- 
ogy, the cores of the interacting clusters must overlap. 
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Similarly, for a given total mass ratio, the relative ram 
pressure experienced by the clusters is determined by 
the relative central gas densities. The central gas den- 
sity ratio {rieoi / neo2) was varied from 0.5 to 2.0. The 
goal was to have enough ram pressure to create the 
desired X-ray core morphology but not so much that 
the subcluster penetrates the primary core. 

The central gas densities are determined by a com- 
bination of the global gas fraction (fb) and the distri- 
bution of the gas relative to the dark matter (i.e., the 
/3-parameter). We do not attempt to determine the 
global gas fraction from the data. Rather, we ar- 
bitrarily choose the more massive cluster to have 
fb=0.12 and /3=0.75 (See §p^). Both values are 
within the acceptable range defined by observation 
(e.g., Jones & Forman 1984; White & Fabian 1995; 
Lubin & Bahcall 1993). We then varied these pa- 
rameters for the less massive cluster, also within the 
observed range, in order to obtain the desired central 
gas density ratio. 

Below we briefly summarize the initial cluster pa- 
rameters. We further discuss the specific effects of 
parameter space in §^ 

4.4. A Brief Model Summary 

Table 2 contains the initial cluster parameters. 
These clusters were allowed to merge under the in- 
fluence of their mutual gravity. They are given an 
initial velocity of 270 km s~^ parallel to the line of 
centers and 100 km s^^ perpendicular. This results 
in a slightly off-axis merger with an impact parame- 
ter of ~120 kpc and a final impact velocity of 2700 
km s""'^. Both of these values are relative to the re- 
spective centers of mass of each cluster. The simula- 
tion most closely resembles A754 at ^0.3 Gyrs after 
closest approach. Although the relative masses of the 
two clusters is poorly constrained, it is important that 
the SE component (the primary cluster in our model) 
have the greater central gas density initially. It is also 
important that the NW component (secondary clus- 
ter in our model) have a relatively more concentrated 
mass distribution. Below, we discuss the justification 
for the choice of these parameters. 

5. A COMPARISON WITH THE DATA 

Figure |^ contains plots of four quantities overlaid 
with the X-ray surface brightness contours. They in- 
clude (a) the projected, emission-weighted ICM tem- 
perature, (b) the dark matter distribution, (c) the gas 



velocity in the plane of the merger, and (d) a dynam- 
ically inert quantity used to trace the two gas com- 
ponents, also within a slice taken in the plane of the 
merger. In each case, the data is taken from the high 
resolution simulation (see but it represents only 
a small fraction of the total volume (80 x 80 zones, 
in projection). The X-ray emissivity was calculated 
using the Mewe-Kaastra-Liedahl emission spectrum 
for optically-thin plasmas supplied with XSPEC (Ar- 
naud 1996). The X-ray surface brightness is simply 
a line-of-sight integration of the volume emissivity. 
Similarly, the emission- weighted, projected ICM tem- 
perature is simply a line-of-sight integration of the 
product of the temperature within a given zone and 
the X-ray emissivity in that same zone divided by the 
total emissivity along a given line-of-sight. The dark 
matter distribution is the projected particle density. 
The merger is assumed to be exactly in the plane 
of the sky. The epoch depicted is 0.3 Gyrs after the 
closest approach of the respective dark matter centers 
of mass. The smaller of the two clusters has moved 
left to right passing below the core of the more mas- 
sive cluster. Hereafter, we refer to the lower left dark 
matter concentration as the primary cluster, and to 
the upper right dark matter concentration as the sub- 
cluster. When comparing to the observational data, 
lower left is SE, upper right is NW. Table 3 shows a 
comparison of some of the global cluster parameters. 

5.1. X-ray Morphology and Hydrodynamics 

The simulated X-ray surface brightness morphol- 
ogy at 0.3 Gyrs is represented by contours in each of 
the four panels in Figure I Using our model, we are 
able to reproduce the basic features of A754's X-ray 
morphology (§2.2., Figure These include, the bar- 
shaped peak in the X-ray emission, the displacement 
of the X-ray peak relative to the galaxies/dark mat- 
ter, and the overall east-west elongation of the X-ray 
emission. The relationship between the orientation of 
the emission peak to the overall X-ray surface bright- 
ness distribution represent an isophotal twisting of 
nearly 90 degrees. We now identify the underlying 
processes responsible for the morphological features 
exhibited by the simulated data. 

The flattening of the X-ray peak on the eastern 
side is a result of the extreme ram pressure {pv'^) of 
the subcluster on the primary. At the time of closest 
approach (as defined by the respective dark matter 
centers of mass) the relative velocity is ^^2500 km 
s^^. Since the merger is slightly off-axis, a portion of 
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the subcluster (moving east to west) essentially runs 
into a wall of gas (i.e., the core of the primary clus- 
ter) with sufhcient ram pressure to flatten the dis- 
tribution and displace it relative to the dark matter. 
Figure ^ shows the basic flow pattern within a 2-D 
slice in the plane of the merger. The maximum gas 
velocities at this epoch are associated with the resid- 
ual infalling subcluster gas located SE of the X-ray 
peak. These gas velocities are seen to approach 1800 
km s""'^ (relative to the primary cluster dark matter 
distribution), but they abruptly decrease in magni- 
tude (<1000 km ) as they pass through a shock 
to the SE of the core after which the infalling gas is 
deflected by the core toward the SW and NE. StiU, 
there is sufflcient ram pressure to displace the X-ray 
peak from the gravitational potential minimum by a 
signiflcant distance. 

ZZ95 flnds the peak of the X-ray emission to be 
displaced by 3.8' from the peak of the SE galaxy 
concentration (assumed to be the minimum of the 
gravitational potential). Looking at the more statis- 
tically robust centroids of these distributions, the off- 
set is somewhat less, ^^3.0' . Our model shows a dark 
matter/X-ray offset of 100 kpc which is ~1.5' assum- 
ing Ho =65 km s~^ Mpc^^. Although somewhat less 
than the observed value, the simulated offset is in the 
same direction as that seen in the data and of compa- 
rable magnitude. A larger offset may be obtained by 
increasing the ram pressure of the secondary relative 
to the primary. However if it is increased too much, 
the subcluster will actually penetrate the core of the 
primary and the bar-like structure will be destroyed. 
Similarly, decreasing the impact parameter may also 
increase the offset. 

Since the merger is slightly offaxis, a portion of the 
subcluster misses the primary core and meets with 
signiflcantly less resistance as it passes to the south 
of the X-ray peak. Figure ^ shows the location of the 
respective gas components relative to the X-ray emis- 
sion. The red indicates subcluster gas while black in- 
dicates primary cluster gas. In this analysis, we have 
used a dynamically inert quantity to trace the indi- 
vidual gas components within the velocity fleld. We 
flnd that the X-ray peak is composed almost entirely 
of gas originating with the primary cluster. The bulk 
of the subcluster gas is located SE of the X-ray peak. 
Very little mixing of the two components has occurred 
at this time. The leading edge of the subcluster gas 
appears as the westward extension of red. As the sub- 
cluster gas passes by the core, it creates a shear which 



drags both primary and subcluster gas into an ex- 
tended tail along the bottom edge of the X-ray peak. 
The overall morphology of the remnant X-ray core is 
seen to have an angled or L-shaped distribution. This 
is reminiscent of a feature identifled in the wavelet 
analysis of A754's X-ray surface brightness by Slezak 
et al. (1994) (Figure 16, therein). Along the interface 
between the two gas components, we note a spray of 
gas both to the NE and SW, essentially perpendicular 
to the merger axis. The velocity of the spray is ^1000 
km s~^ near the X-ray peak but is seen to broaden 
and decelerate as it exits the core. 

Comparison of the location of the subcluster gas 
component in Figure ^ with the location of the dark 
matter distribution in Figure |^b reveals that the sub- 
cluster gravitational potential (NW) has been com- 
pletely stripped of its original gas content. A similar 
stripping was noted by Pearce et al. (1994). This 
is not to say that the remnant subcluster is devoid 
of gas. Rather, it drags along gas originally associ- 
ated with the primary causing it to interact strongly 
with gas infalling from the outer regions of the clus- 
ter. The interface between the expanding bow shock 
and infalling gas is most noticeable in the temper- 
ature distribution (Figure ||a, § |5.2.[ ) and in the gas 
velocity vectors (Figure ^) where a sharp transition 
is seen in the NW between gas moving to the NW 
at up to 1000 km s^^ and gas moving to the SE at 
nearly 700 km s^^. Regarding the X-ray morphology, 
it is the gas being dragged outward by the subclus- 
ter potential that accounts for the overall (SE to NW) 
elongation of the X-ray surface brightness. The impli- 
cations for gas temperature morphology are discussed 
below (jj5?2l). 

Of course the gasdynamics in Figure |^c would not 
be directly observable since they are in the plane of 
the sky. In Figure ^ we show the line-of-sight gas ve- 
locities (vertical axis) in an east-west scan (horizon- 
tal axis) which intersects the X-ray surface brightness 
peak. That is. Figure || represents gas velocities in a 
plane perpendicular to the plane of the sky. Evident 
in this figure is a general expansion of the cluster gas. 
The largest velocities (~1250 km s~^ ) are in the re- 
gion of the radial spray described above. We also see 
expansion velocities near 500 km s~^ in the region 
of the bow shock (right hand side). These velocities 
should be resolved by ASTRO-E provided the emis- 
sivity is great enough in the high velocity regions. Of 
course in an actual spectroscopic observation, the dy- 
namical structure seen in Figure would be projected 
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along the line-of-sight and emission- weighted. 

5.2. Temperature Distribution 

Figure |^ shows the projected, emission-weighted 
temperature distribution (color) upon which is su- 
perimposed the X-ray surface brightness distribution 
(contours). The basic temperature morphology con- 
sists of two hot spots. The smaller, lower tempera- 
ture hot spot located somewhat south and east of the 
X-ray peak is formed when residual gas from the sub- 
cluster collides with gas stalled near the X-ray peak. 
As such, this feature is not obviously apparent in the 
observational data although it would be more consis- 
tent with the data if it were a somewhat lower tem- 
perature and/or further S/SW. Our parameter space 
survey shows it to be a strong function of both the 
relative gas content of the clusters and the impact 
parameter. Decreasing the gas content of the sub- 
cluster relative to the primary decreases the hot spot 
temperature relative to the surrounding gas by de- 
creasing the ram pressure and limiting the residual 
infall. Decreasing the impact parameter moves the 
hot spot to the north and increases its temperature 
relative to the surrounding gas. Therefore, a better 
fit to the data may be possible by slightly increasing 
the impact parameter and decreasing the gas content 
of the subcluster. However, this must be balanced 
by other features. Specifically, we need sufficient gas 
in the subcluster to create the ram-pressure flatten- 
ing and displacement of the X-ray peak (§5.1.). This 



would argue for a more concentrated gas distribution 
in the subcluster. 

The second hot spot appears as a large slightly 
asymmetric arc near the western edge of the X-ray 
emission. Figure The peak projected, emission- 
weighted temperature is greater than 19 keV and 
largely coincides with the subcluster dark matter dis- 
tribution. As the subcluster impinges on the pri- 
mary a bow shock forms which produces a conical 
sheet of extremely hot gas. As the merger progresses, 
the shock heated gas is driven through the primary. 
When viewed in projection, it appears as the arc that 
we see here. Although the subcluster is stripped of its 
initial gas content, the dark matter potential remains 
intact and drags a significant amount of primary clus- 
ter gas with it. This gas continues to interact with 
gas in the outer regions of the primary and helps to 
maintain the temperature of the original shock heated 
gas. Coincident with the arc of hot gas is a disconti- 
nuity in the gas velocities which approaches 1700 km 



Another interesting temperature feature is the band 
of cooler gas (^7.5 keV) that runs along the X-ray 
peak. This is a combination of pre-shocked primary 
cluster gas (initially 6.7 keV) and subcluster gas that 
is deflected by the core and forced outward in the ra- 
dial spray mentioned above (see §5.1,). As the spray 



expands and decelerates, the temperature drops to 
~6.5 keV at 0.5 Mpc to the NE. 

Figure ^ shows a quantitative point-by-point com- 
parison of our model with the ASCA-based temper- 
ature map of A754 (HM96). Since we find a general 
consistency between the ASCA-based and ROSAT- 
based (HB95) temperature maps, we will only com- 
pare directly with the former. For the purposes of this 
analysis, we have placed a grid similar in scale and 
location (with respect to the X-ray emission) to the 
one used by HM96. Within each numbered region, 
we have calculated the mean projected, emission- 
weighted temperature in the model. We then compare 
these to the values observed in the corresponding re- 
gions. Error bars on the observed values are 90% con- 
fidence levels. Error bars on the model represent 1-a 
about the mean. Data points have been slightly offset 
to minimize confusion. In each instance, the model 
is statistically consistent with the observed temper- 
atures. The largest absolute discrepancy is found in 
region 7. We do find gas at ~19.5 keV in this region, 
however because of the exact location of the grid, it 
becomes somewhat diluted with adjacent, cooler gas. 

In these data, we recognize three basic features. 1) 
The negative gradient from region 1 to region 3. 2) 
The relatively isothermal regions 4,5, and 6. 3) The 
negative gradient from region 7 to region 9. Each of 
these features is well-represented in the model. The 
largest discrepancy is found in regions 4-6 where the 
model shows a slightly positive gradient and the data 
shows a somewhat stronger negative gradient because 
of the relatively hot region 4. Still, they are statis- 
tically consistent. As was suggested above, better 
agreement in region 4 might be obtained if the hot 
spot near the X-ray peak (Figure |^) were further S- 
SW. 

Of course, the X-ray telescope only allows us to 
see the projected, emission-weighted temperature. In 
doing so, much of the detailed structure of and ex- 
tremes in temperature can be hidden. Figure |^ shows 
the zone-by-zone distribution of the true (i.e., non- 
emission weighted) model gas temperatures within a 
cube comparable in dimension to that delineated by 
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the 9 regions in Figure |[ Note the wide range in 
temperatures (5-22 keV). Although the bulk of the 
volume is dominated by gas less than 10 keV, the dis- 
tribution is strongly skewed toward higher tempera- 
tures. Fully, 29% of the volume is in gas greater than 
11 keV. The median temperature is 9.9 keV, the mean 
temperature is 10.9 keV. Compare these values to the 
mean emission- weighted temperature of 9 keV and the 
initial primary cluster temperature of 6.7 keV. 

6. FUTURE EVOLUTION 

Figures ^ and ^ depict the model evolution at 2.75 
and 5.25 Gyr after closest approach. These are the 
same quantities represented in Figure]^ (T=0.3 Gyr) 
and can therefore be compared directly. Considerable 
evolution is evident. The subcluster dark matter dis- 
tribution has proven to be rather resilient. Not only 
does it survive the initial passage, but a couple of 
others as well. Through examination of the particle 
group velocities, we identify subsequent merger events 
at 1.6 and 2.5 Gyr. Although far less dramatic than 
the initial merger, these events are still significant in 
that they continually stir the ICM preventing a rapid 
relaxation. At 2.75 Gyrs, the system is experiencing 
a maximum separation which results in the peanut- 
shaped morphology of the dark matter distribution 
seen in Figure ^ (color). It should be noted that 
survival of substructure in the dark matter will de- 
pend on the relative masses of the subclusters as well 
as their relative concentration. As the mass ratio de- 
creases (approaching 1:1), the subcluster will survive 
longer. Also, chances of survival will be enhanced as 
the subcluster mass distribution becomes more con- 
centrated and as the impact parameter increases. 

The continued dynamical evolution is still appar- 
ent in the X-ray surface brightness morphology (all 4 
panels, Figure ^). Although considerably more sym- 
metric (i.e., relaxed) than during the epoch we find 
most consistent with A754, there is still a significant 
extension to the SE which is abruptly flattened. Com- 
parison with the velocity vectors in Figure ||c, indi- 
cates that the extension results from a southerly flow 
of ~650 km s~^. The flattening occurs in a region 
of strong shear where the southerly flow is abruptly 
cutoff by a WSW flow of -700 km s'^. In fact, the 
WSW flow is part of a more global circulation that is 
seen to extend unbroken (at least in the plane of the 
merger) around the entire cluster. In this sense, the 
region of shear described above is the location where 



the flow meets itself. One by-product of the global cir- 
culation is the mixing of the subcluster and primary 
cluster gas components. Figure ^ demonstrates how 
the global circulation has drawn the subcluster gas 
(red) around the remnant core. We examine the rate 
of gas mixing more quantitatively in §6.3. Finally, the 
temperature extremes seen in A754 have been reduced 
significantly by 2.75 Gyrs (Figure ||a). At this time, 
the cluster would be considered largely isothermal, at 
least within the central region. To the NW, we find 
large regions of relatively cool gas (<6 keV) which are 
likely remnants from the initial merger. Although ra- 
diative cooling was not included in these simulations, 
we do find regions where the gas cools via adiabatic 
expansion such as behind the initial bow shock. 

Between 2.75 and 5.25 Gyr, there continues to be 
minor "sloshing" within the gravitational potential. 
The relative dark matter center of mass velocities 
peak at only 100 km s^^. At 5.25 Gyr, the X-ray 
surface brightness morphology is quite regular with a 
slight elongation from SE to the NW (Figure 0). A 
similarly relaxed dark matter distribution is apparent 
in Figure l^b (color). Here, the dark matter distribu- 
tion appears somewhat more elongated than the gas 
distribution. Like the gas distribution, the dark mat- 
ter major axis is oriented SE to NW reflecting the axis 
of the secondary mergers. Unlike the gas distribution, 
the dark matter can shed angular momentum through 
free streaming particles. Consequently, the secondary 
mergers discussed above are largely head-on. 

Although temperature inhomogeneities persist, the 
overall temperature distribution is noticeably more 
regular (Figure ^ than at earlier epochs. The cool 
regions noted at 2.75 Gyr have vanished while a radial 
temperature gradient peaking at 11 keV in the core 
has developed. Pearce et al. (1994) also note a cen- 
trally peaked temperature distribution. The relative 
heating in the core results from the continual stirring 
of the relatively high density gas in the core by fluc- 
tuations in the gravitational potential minimum as 
particles continue to stream back and forth along the 
merger axis. In the absence of radiative cooling, there 
is no way to dissipate the heat. We also find that the 
region NW of the cluster core is somewhat warmer 
(—1-2 keV) than the corresponding region to the SE. 
The overall temperature morphology at this time is 
similar to that of Triangulum Australis (Markevitch, 
Sarazin & Irwin 1996). 

Largescale gasdynamics (both bulk flows and rota- 
tion) persist even as late as 5.25 Gyr. Still evident by 
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the velocity vectors in Figure 0c is a general counter- 
clockwise circulation in the merger plane. Comparing 
to Figure |6|3, we note that the flow is less coherent and 
of significantly reduced magnitude. The highest ve- 
locity gas ('~850 km s~^ ) is located in the extreme 
SE corner of Figure |^. Within the central cluster, 
gas velocities range up to 600 km s~^, but the aver- 
age in nearer 300 km s"'^. As mentioned above, the 
gas is less efficient than the dark matter at shedding 
angular momentum. Over time, angular momentum 
is transferred via mixing to adjacent gas components. 
Similarly, kinetic energy is transferred from one com- 
ponent to the next and ultimately converted into heat. 
The relative inefficiency of this process is similarly ap- 
parent in the slow rate of mixing between primary and 
subcluster gas components, particularly in the outer 
regions of the merger remnant (Fig ure We ad- 
dress mixing in more detail below, §|6.3.|. 



6.1. Evolution of Dynamical Pressure Sup- 
port 

Since dynamical information regarding the X-ray 
emitting gas is non-existent at this time, the analy- 
sis of clusters has always assumed them to be ther- 
mally supported (i.e., in hydrostatic equilibrium). 
Several recent cluster analyses have argued for ad- 
ditional pressure components (e.g., dynamical, mag- 
netic fields) in order to explain the apparent discrep- 
ancies between X-ray, virial and gravitational lensing 
based cluster mass estimates (e.g., Loeb & Mao 1994, 
among others) . Using our model of A754, we can de- 
termine the relative significance of dynamical pressure 
support throughout the merger evolution. 

Figure || shows the evolution of Pk / Pt {Pk = pv^ 
and Pt = kT) within spherical volumes (r=0.4 and 
1.5 Mpc) centered on the gravitational potential min- 
imum. Initially, the dynamical pressure is negligible 
although it grows steadly as the subcluster impinges 
on the primary. A maximum is reached shortly before 
closest approach (T=0) when the dynamical pressure 
is greater than 60% of the thermal pressure within 
0.4 Mpc. Secondary maxima occurring near 2 and 3 
Gyr correspond to subsequent passages by the sub- 
cluster dark matter remnant. As might be expected, 
the larger volume (r=1.5 Mpc) experiences fewer ex- 
tremes than the core which may explain the relative 
heating in the core. After the initial impact, the dy- 
namical pressure is typically '^20% of the thermal 
pressure but decays to 5-10% after 4.5 Gyrs. There- 
fore, late in the merger evolution, gasdynamics play 



a limited role in supporting the expanded gas core. 
This is consistent with Pearce et al. (1994). 

Figure || has implications for X-ray based cluster 
mass estimates. Neglecting the dynamical component 
of the pressure support will cause an underestimation 
of the cluster mass by a fraction 1-{1+{Pk / Pt))~^ ■ 
Based on this analysis, we would predict errors of 
~20-40% or more during the 2 Gyrs immediately af- 
ter a merger and errors of 10-20% at later times. Note 
also the radial dependence. Mass estimates based on 
only the central regions of the cluster will be more 
severely affected. These results are consistent, at least 
in magnitude, with previous studies (e.g., Evrard, 
Metzler & Navarro 1996; Schindler 1996; Roettiger 
et al. 1996) although none of these studies showed 
evidence for systematic underestimation of the clus- 
ter masses indicating that other effects (e.g., projec- 
tion effects) and the details of the observational anal- 
ysis are also important. By underestimating the to- 
tal mass, the baryon fraction would be overestimated 
by the same factor. Although significant, the errors 
noted here are not enough to account for the apparent 
inconsistency between measurements of the baryon 
fraction in clusters and that derived from primordial 
nucleosynthesis models that assume $1=1 (White et 
al. 1993). Of course lowering the assumed value of 
also reduces the perceived discrepancy. 

6.2. Angular Momentum Evolution 

The offaxis merger imparts considerable and sus- 
tained angular momentum upon the gas distribution 
in the merger remnant (Fig|^). Although the gas dis- 
tribution is not rotationally supported, it is slightly 
less prolate than the dark matter distribution where 
angular momentum is dissipated more readily by free- 
streaming particles. Figure |^ shows the evolution of 
total angular momentum density (Ltot) within con- 
centric spheres of radius 0.4 and 1.5 Mpc centered 
on the gravitational potential minimum. We define 

Ltot = [S?(i.)']°''/^. and = Sf (r, x p,v,)i 
where Vj is the distance from the potential minimum 
to the jth zone, pj and Vj are the gas density and ve- 
locity, respectively, within the jth zone, and Nz is the 
number of zones within the volume. The sum over i 
represents the sum over the Cartesian coordinates. 

The angular momentum evolution in Figure]^ closely 
mimics the evolution of the non-thermal pressure sup- 
port (see §3.1., Figure @). The abrupt rise at T=0 (R 



< 0.4 Mpc) occurs as the subcluster passes near the 
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primary core. The subsequent decay occurs as the 
subcluster leaves the sphere only to return at 1.5-2.0 
Gyr. Beginning at ~3.5 Gyr, the angular momentum 
within 0.4 Mpc is seen to decay steadily as the core 
relaxes and angular momentum is transferred via mix- 
ing to larger radii. Within the larger volume, the total 
angular momentum remains fairly constant, peaking 
slightly at 1.5-2.0 Gyr. Once again, this shows that 
although the canonical core sound crossing time (and, 
therefore its relaxation time) is of order 1 Gyr, true 
relaxation of the merger remnant does not really be- 
gin until after ~3 Gyr. The longer relaxation time 
is important to the interpretation of substructure fre- 
quency. A longer relaxation time allows for greater 
consistency between a high substructure frequency 
and a low fl, the density of the universe relative to 
the critical density (Boute & Xu 1997). 

6.3. Mixing 

Using a dynamically inert quantity evolved within 
the velocity field, we have traced the rate at which the 
primary and subcluster gas components mix. In addi- 
tion to the non-isothermal temperature distributions 
observed in some clusters, Arnaud et al. (1994) have 
also noted an inhomogeneous abundance distribution 
in Perseus. Since considerable scatter is observed in 
global cluster abundances (Mushotzky & Loewenstein 
1997), it may be reasonable to expect inhomogeneous 
distributions within individual clusters to arise as a 
result of mergers between clusters of differing initial 
abundances. It then becomes important to quantita- 
tively assess the rate at which the remnant gas dis- 
tribution becomes homogeneous. A cursory examina- 
tion of Figures |^, and |^ indicates that mixing 
is slow in any case but even more so at larger radii 
where the gasdynamics are less extreme. Even after 
5.25 Gyrs, there is a noticeable patchiness in the gas 
distribution. This result may be sensitive to cosmo- 
logical infall which we do not simulate. 

Figure |l^ presents the degree of mixing within 
spheres of radius 0.4 and 1.5 Mpc. Plotted is the frac- 
tion of zones within the corresponding volumes that 
are at least 10% mixed. Mixing does not begin un- 
til after the time of closest approach (T=0) at which 
point the degree of mixing within both 0.4 and 1.5 
Mpc increases steadily with time. The level of mix- 
ing is consistently lower in the larger volume although 
the rate of mixing within the two volumes track well 
until about 3-4 Gyrs after the initial passage. At this 
time, mixing within the 1.5 Mpc sphere effectively 



stops while the rate of mixing within 0.4 Mpc actu- 
ally increases. Recall from §3.2. that it is at 3.5 Gyr 
when the ordered circulation in the core gives way to 
random motions thus increasing the mixing efficiency. 
After 5 Gyrs, 90% of the zones within 0.4 Mpc are at 
least 10% mixed which compares to only 45% within 
the larger volume. This analysis would seem to in- 
dicate 1) a tendency to find more homogeneous gas 
distributions within the cores of clusters, 2) the po- 
tential for abundance gradients within clusters having 
evolved sufficiently since their last significant merger, 
and 3) the potential for patchiness in the distribution 
of metals in the outer regions of clusters. Of course a 
patchiness may also result from a local enhancement 
of metals by supernovae in individual galaxies. 

7. SUMMARY 

We have presented a dynamical 3-dimensional model 
of A754 which can explain many of its observed mor- 
phological peculiarities. In our model, A754 is the 
result of a very recent (<0.3 Gyr), nearly in the plane- 
of-the-sky merger between two clusters having a total 
mass ratio that is likely less than 2.5:1. The merger is 
slightly off-axis, with an impact parameter of 120 kpc 
or about half the initial primary cluster core radius. 
The final impact velocity is ^2500 km s^^. Although 
the model is poorly constrained by the limited ob- 
servational data, our assumed initial conditions for 
the merger are consistent with that data. The fact 
that the observed characteristics of A754 can be re- 
produced so well suggests the actual properties of the 
merging clusters in A754 cannot be too different from 
those adopted here. Moreover, we believe that the 
dynamics described here are largely representative of 
slightly off-axis mergers. 

In general, the results of these simulations are con- 
sistent with other studies of major cluster mergers 
( Roettiger et al. 1993, 1996, 1997a; Schindler & 
Miiller 1993; Pearce et al. 1994). In particular, they 
reveal how gravitational energy is transferred to the 
ICM and dissipated via shocks, resulting in a heat- 
ing and expansion of the ICM relative to that of the 
dark matter, and high velocity bulk flows within the 
gas. It should be noted however that in this study, 
as well as many others, potentially important physical 
processes (e.g., radiative cooling, thermal conduction, 
magnetic fields, etc.) have been ignored. Moreover, 
we study the merger of isolated subclusters in hy- 
drostatic equilibrium rather than mergers within an 
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evolving largescale structure. On the other hand, this 
allows us to focus our numerical resources on the dy- 
namics of the merger itself. In addition, we are able 
to systematically identify and analyze the physical 
processes which govern the dynamics during mergers. 
Thus, our results have important implications for the 
dynamics of off-axis mergers in general. 

More specifically, the main conclusions of this 
study include: 

1) Major merger events can account for many if 
not all of the morphological irregularities observed in 
A754 as well as other well-resolved clusters. The ex- 
treme temperature inhomogeneities of the type ob- 
served in A754 are currently the most direct indica- 
tor of recent dynamical evolution and are particularly 
useful for ruling out more subtle projection effects. 
The simple projection of two clusters along a com- 
mon line-of-sight cannot account for the sharp tem- 
perature transitions, due to shocks, nor the extreme 
nature of the observed temperatures (19 keV). 

2) Relaxation of the merger remnant is delayed 
beyond the canonical sound crossing time by subse- 
quent passages of the subcluster's remnant dark mat- 
ter. The extended relaxation time allows for greater 
consistency between a high frequency of substructure 
in clusters of galaxies and a low fl universe (Buote & 
Xu 1997). 

3) Major merger events can result in significant 
non-thermal pressure support within the remnant gas 
distribution. The existence of a non-thermal pressure 
component may seriously affect X-ray based cluster 
mass estimates which often employ the assumption of 
hydrostatic equilibrium. Errors are lower when the 
mass estimates are made at larger radii, but they can 
still approach 30% within 1.5 Mpc during the early 
stages of a merger. At later times, errors are more 
typically 15%. These results are consistent in magni- 
tude with Evrard et al. (1996) and Schindler (1996). 

4) Off-axis mergers can impart significant angu- 
lar momentum on to the remnant gas distribution re- 
sulting in a sustained rotation about the cluster core. 
The magnitude of rotation (and bulk flows) seen in 
this model should be observable with ASTRO-E. Al- 
though much of the gasdynamics in A754 merger are 
not directly observable since they occur in the plane 
of the sky, we do predict that a general expansion 
of the cluster ICM may be observable provided the 
emissivity of the high velocity gas is great enough, see 
§5.1. Since the time period between major mergers is 



likely comparable to or less than the observed relax- 
ation time scales, the gasdynamics within the cluster 
may largely reflect those imparted by the last major 
merger. 

5) Mixing of premerger gas components is an inef- 
ficient process and may explain the observed patch- 
iness and gradients in both temperature and abun- 
dance distributions (e.g., Perseus, Arnaud et al. 1994; 
Coma, Honda et al. 1996). Mixing is most efficient 
in the core of the cluster where the gasdynamics are 
most significant. 
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Tal)l(> 1. ()l)S(>r\"(Hl Sul)c-lusl(>r Gala.x\- Paraiiiclt^rs 





Vr (km 




(km 










SE^ 




SE 


NW 


SE 


NW 


F86* 


16145±160 


16165±217 




839±^^i 


21 


15 


EM925 


16387 


16007 


698 


521 


12 


10 




16320±132 


16432±191 


830±i?^ 


970±lil 


41 


27 



Table 1: ^ Number of galaxies associated with the subcluster ^ Southeastern subcluster, ^ Northwestern subcluster, 
* Fabricant et al. 1986, ^ Escalera & Mazure 1992, *^ Zabludoff & Zaritsky 1995 



Tal)lc 2. Initial C'lusli^r Paraiutncrs 



Cluster^ 
ID 


(10" Mq) (keV) 


(km s-i) 


rl 
(kp"c) 


Jg 




"-eo 

(10-3cm-3) 


(km s-i) 


(kpc) 


1 


8.0 6.7 


785 


220 


0.12 


0.75 


1.55 


2500 


120 


2 


3.2 3.1 


526 


135 


0.06 


0.72 


1.01 







^ Cluster ID. ^ Total Mass R<3 Mpc. ^ Temperature. ^ Velocity Dispersion. 

^ Core Radius. ^ Global Gas Fraction, by mass. ^ /3-paramctcr. 

^ Central Gas Density. ^ Impact Velocity (dark matter). ^° Impact Parameter. 



Table 3. Post-Merger Comparison 





N,(SE)^ 


Ng(NW)^ 


a^(SE)^ 
(km s-i) 


a^(NW)4 
(km s-i) 


< Te >^ M«,, 
(keV) (10" Mq ) 


Observed 


41 


27 


830±i5'^ 


970±|§t 


9.0 10.48±3.27 


Simulated 


50000 


20000 


1034±105 


830±145 


9.1 9.5 



^ Galaxy/particle count in SE subcluster (ZZ95). 
^ Galaxy/particle count in NW subcluster (ZZ95). 
3 Velocity dispersion in SE (ZZ95). 
^ Velocity dispersion in NW (ZZ95). 
^ Mean temperature (HM96). 

^ Total mass within -^2 Mpc (Escalera et al. 1994, virial mass). 
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FIGURE CAPTIONS 



Fig. ^ ROSAT X-ray surface brightness (solid contours) and projected, galaxy density (dashed contours). The 
galaxy distribution is derived from Fabricant et al. (1986) and Dressier & Shectman (1988). In all, 142 galaxy 
positions were smoothed with 2.5' Gaussian. The X-ray surface brightness is smoothed with a 40" Gaussian. 
Contours are at 0.08, 0.16, 0.24, 0.32, 0.4, 0.48, 0.56, 0.64, 0.72, 0.80, 0.90, 0.99 of the peak. 

Fig. ^ Simulated data at 0.3 Gyr after closest approach. Each panel is 2 Mpc x 2 Mpc. The smaller of the 
two clusters (subcluster) entered from the left and is moving to the right, a) X-ray surface brightness (contours) 
and projected, emission- weighted temperature (color). Note the color bar above the panel, b) Projected dark 
matter particle distribution (color) c) Gas velocity within the plane of the merger. Vectors are spaced at 75 kpc 
( one third the resolution of the simulation) and scaled to the maximum velocity indicated below the panel, d) 
A dynamically inert quantity used to trace the subcluster gas component (red) as it mixes with the primary gas 
(black). Intermediate colors indicate the degree of mixing as shown in the color bar below the panel. The X-ray 
contours are at 0.03, 0.08, 0.16, 0.24, 0.32, 0.4, 0.48, 0.56, 0.64, 0.72, 0.80, 0.90, 0.99 of the peak. These are the 
same as in Fig. ^ with the exception of the lowest contour. 

Fig. ^: Line-of-sight gas velocities in an east-west scan that intersects the X-ray surface brightness peak. The 
vertical axis represents the observer's line-of-sight looking into the cluster. The horizontal axis represents an east- 
west scan across the projected face of the cluster. Contours are ±250, ±500, ±750, ±1000, ±1250 km s~^. Negative 
velocities (approaching) are represented by dashed contours. Positive velocities (receding) are represented by solid 
contours. Tick marks represent 0.125 Mpc. The plot dimensions are the same as the panels in Fig. ^. In a 
spectroscopic observation of the cluster, the velocity structure seen in this figure would be emission-weighted and 
projected onto the plane of the sky. Consequently, much of the detailed structure may be unobservable. 

Fig. ^ A comparison with the HM96 temperature map. a) Simulated X-ray surface brightness (contours) and 
grid similar in scale and location (relative to the X-ray image) as that used by IIM96. b) Simulated projected, 
emission-weighted temperatures (*) within the numbered regions in (a) compared directly with the corresponding 
temperatures observed by IIM96 (o). 

Fig. |: A histo gram showing the zone-by-zone distribution of temperatures within the volume delineated by 
the nine regions in Figure ^. 

Fig. ^ Same as Figure |2| only at 2.75 Gyr after closest approach. 

Fig. ^ Same as Figure |2| only at 5.25 Gyr after closest approach. 

Fig. ^: The evolution of the ratio of dynamical to thermal pressure support within spheres of radius 0.4 Mpc 
(solid line) and 1.5 Mpc (dashed line) centered on the gravitational potential minimum. Times are relative to the 
time of closest approach. This analysis was performed on the low resolution simulation (§^. 

Fig. ^: The evolution of the angular momentum density per zone relative to the gravitational potential minimum 
within spheres of radius 0.4 Mpc (solid line) and 1.5 Mpc (dashed line) centered on the potential minimum. Times 
are relative to the time of closest approach. This analysis was performed on the low resolution simulation ([js]). 

Fig. 0: The fraction of zones within spheres of radius 0.4 Mpc (solid line) and 1.5 Mpc (dashed line) centered 
on the gravitational potential minimum which exhibit at least 10% mixing. That is, the gas within that zone 
consists of at most 90% of either primary or subcluster gas. This analysis was performed on the low resolution 
simulation ([|3]). 
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